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Effects of anisotropy on the geodesic acoustic mode (GAM) is analyzed by using gyro-kinetic 
equations applicable to low-frequency microinstabilities in a toroidally rotating tokamak plasma. 
Dispersion relation in the presence of arbitrary Mach number M, anisotropy strength cr, and the 
temperature ration r is analytically derived. It is shown that when cr is less than 3 + 2r, the 
increased electron temperature with fixed ion parallel temperature increases the normalized GAM 
frequency. When cr is larger than 3-|-2r, the Increasing of electron temperature decreases the GAM 
frequency. The anisotropy a always tends to enlarge the GAM frequency. The Landau damping 
rate is dramatically decreased by the increasing r or cr. 
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I. INTRODUCTION 

Tokamak plasmas are usually described by isotropic 
equilibria, even in the low-collisional or collisionless 
limiti . It is a reasonable approximation for Ohmically 
heated discharges but becomes quite unrealistic in the 
presence of intense auxiliary heating. It is known that the 
neutral beam injection (NBI) heating, ion cyclotron reso¬ 
nance heating (ICRH), and electron cyclotron resonance 
heating (ECRH) can produce strong plasma anisotropy^. 
More specifically, NBI generally creates anisotropic equi¬ 
libria with p|| > p_L, while the ICRH typically drives 
perpendicular dominant anisotropy with p± > p||, where 
P\\{p±) is the parallel (perpendicular) pressure with re¬ 
spect to the magnetic field^. In some experiments at 
JETl, an ICRH-heated minority ion population with 
strong temperature anisotropy was observed. In MAST, 
the beam pressure ratio p±/p\\ — 1.7 was achieved dur¬ 
ing the NBI heating^. Plasma equilibrium with strong 
anisotropic effects was also found in the low density dis¬ 
charges of the Large Helical Device (LHD) with powerful 
tangential NBI^. 

The theory of tokamak anisotropic equilibria and 
stability has been approached long ago^“— by many 
authorsi^^— devoted to the equilibria by using magneto¬ 
hydrodynamic (MHD) model or the guiding-center the¬ 
ory (see, for example. Refs. 00110), and to the stabil¬ 
ity by using energy principl o^°i^^i^^ . The importance of 
the pressure anisotropy effect on the plasma equilibrium 
in tokamaks was emphasized in Ref. [2^ by implementing 
anisotropic pressure model into the equilibrium recon¬ 
struction code EFIT and applying the algorithm to se¬ 
lected high-performance discharges of the tokamaks JET 
and Tore Supra. It was shown that due to the assump¬ 
tion of isotropic pressure, earlier attempts using the pres¬ 
sure from transport analysis codes often gave unrealis¬ 
tic results^i^i. While the anisotropic pressure model can 
yield a consistent equilibrium, showing a large difference 
in parallel and perpendicular pressures during additional 
heating^^. One basic outcome of anisotropic equilibria is 
that the outward centrifugal shifts of the magnetic axis 


and of the mass density are proved to be increased by the 
anisotropy with p± > p\\ or decreased by p± < p\\ 

The other is that the two pressures pj_ and p|| and plasma 
density p are no longer flux functions and depend on the 
poloidal angle through the magnetic field strength B. 
The importance of the poloidal asymmetry in the dis¬ 
persion relation of geodesic acoustic mode (GAM) was 
checked in Ref. [13, which is limited to the case of zero 
electron temperature. However, in normal conditions, 
there is Te/Ti ^ 1 in the core and Te/Ti ~ 1 at the edge 
of the plasma column, where Tg and Ti are the electron 
and ion temperatures, respectively. It is of great impor¬ 
tance to take into account electron temperature effect. 
Besides, considering that plasma rotation is also driven 
by auxiliary heating so that a rigid anisotropic equilib¬ 
rium should contain the rotation, which motivated the 
present work devoted to the kinetic GAM in an isotropic 
toroidally rotating tokamak. 


GAM is basically a sort of electrostatic perturbation 
with toroidally symmetrical and poloidally nearly sym¬ 
metrical structure^^. It has attracted intense attention 
since theoretically predicted by Winsor et al. in 19682^. 
The physical interpretation of GAM is that the radial 
drift of the perturbed particles distribution due to the 
geodesic curvature effect compensates the polarization 
drift and an oscillation with a frequency around the ion 
transit frequency is excited. By using drift kinetic equa¬ 
tion, the GAM in the plasma with bi-Maxwellian dis¬ 
tribution for ions was investigated and the previous ki¬ 
netic result was recovered when zeroing the anisotropy^. 
However, the equilibrium distribution function used in 
Ref. 0 actually violated the MHD equilibrium condi¬ 
tion. In the present work, we will derive the dispersion 
relation of GAM in an anisotropic tokamak plasma by us¬ 
ing the modified gyro-kinetic (GK) equation applicable 
to low-frequency microinstabilities in a rotating axisym- 
metric plasma and check the effects of ion anisotropy and 
the ratio between electron and ion temperatures on the 
GAM frequency and Landau damping rate. 
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II. EQUILIBRIUM 

Let us consider a large-aspect-ratio tokamak plasma 
with a toroidally symmetric magnetic field B = 

VC X V'f/’, and work in the {r,9,() coordinate system, 
where is the magnetic flux, and C and 0 are the 

toroidal and poloidal angles, respectively. Circular cross- 
section is assumed in our calculation with R = Rq + 
rcosd. The subscript 0 denotes the equilibrium profile 
and the prefix 6 denotes the perturbed one, but the sub¬ 
script is omitted and the equilibrium magnetic field is 
referred to by B directly. Only toroidal rotation is taken 
into account with uq = uit{'>P)R^'S/(^. w = if — uq is the 
particle velocity in the local reference frame moving with 
Mo relative to the lab frame. In the local frame, a stan¬ 
dard anisotropic Maxwellian distribution function with 
two temperatures is 


^0 = ^^(^,9) 


{ 27 rIrrij) 

rpl/2rp 

J|| 




( 1 ) 


is the greatest^^ and the time scale for electron relaxation 
can be much shorter than the ion time scaled, which re¬ 
quires r ^ mf.jmi. As a result, the ion version of Eq. 
® yields the ions equilibrium density as 


III 


-^+- 
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(7) 


As a naturally result of t ^ me/rui, the mass in the 
electron version of Eq. ([3]) is negligible. One then obtains 
the electron number density 

e^n 

ng = 7Ve(^/>)eT^. (8) 


Here, we write $0 = (‘I’o) + ‘I’o, where (d>o) is the mag¬ 
netic surface average and $0 is the part depending on the 
poloidal angle. Using charge neutrality rig = Mq = ng and 
(rig) = (rig), we arrives at 


6^*0 = -T-— In 
1 -I- r 



1 -I- r 




(9) 


where rig is the number density and vrij is the mass of 
species j {j = i,e for ions and electrons, respectively). 
The zeroth-order Fokker-Planck equation reads^ 
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in which qj{= ±e) is the charge of species j. After sub¬ 
stitution of o, the equation above contains terms with 
three types of velocity dependence, proportional to 
m;|, and w^\w']_. In order that the equation be true for 
all w, it is necessary that the coefficients of each of these 
three types of terms vanish separately: 


6 • (Vlnri^ — Vlnrj_) H - b ■ V(gj$Q -|- m^-Mg • Vug) = 0, 

HI 

(3) 

6-vr||=o, (4) 

(5) 


T± 


&-VlnTi+ (l-^)(V-6)=0. 


Eq. (lU indicates ry = T||(^/i), implies that parallel tem¬ 
perature gradient is forbidden in the local thermal equi¬ 
librium stated. According to Eq. ®, the perpendicular 
temperature has the following form 




B-Boi'ipy 


( 6 ) 


Meanwhile, electrons are supposed to have a standard 
Maxwellian distribution function with temperature Tg. 
It is reasonable since the temperature relaxation rate be¬ 
tween electron perpendicular and parallel temperatures 


where r stands for the temperature ratio Te/T\\ and M 
is the Mach number dehned as utR/v^w with the ion 
parallel thermal velocity vt\\ = \/2T\^lm,i. 


III. DISPERSION RELATION 


The perturbed distribution function is determined by 
the modified GK equation applicable to low-frequency 
microinstabilities in a toroidally rotating tokamak, which 

showsii^i^i 

rivi nr) 

SF, = + [1 - JoikrPj)]^^ScP + JoSh,, 

^ ( 10 ) 

in which 6hj is governed by 

d , - ^ X ^ 

— -I- (w\\b + Uo + VD)-y 

= - ^ + “0 ■ X VS^ ■ VF^ 

■’ dE \dt J rujujl 
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+ Jo^b X y5(j) ■ [(M;||6-k Mg) • VMg -k Vug • iw\\b + Uo)]^^. 
OJc 

( 11 ) 

Only perturbed electrostatic potential is taken into ac¬ 
count, which is justified for electrostatic GAM in a low-/3 
plasma. Here, kr is the radial wave number of GAM, 5(j) 
is the perturbed electrostatic potential, Jg is the zeroth- 
order Bessel function, pj = w±/ijjI is the Larmor radius, 
uj^ = qjB/rrij is the gyro frequency, and md is the leading 
order drift velocity. In the lab frame, electric field can be 

expanded as E = E_i -k Eg -I -, where E_i = —Mg x B 

and Eg is related to ‘I’g, displayed in Eq. (I3|). Hence in 
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the local reference frame, particles only feel the potential 
<i>o- As a result, the drift velocity can be expressed as^ 

vd = [(wy + w5_/2)/w^]& X VIni? 

H-rx -^V3>o + fto ■ Vmo + u>|| (6 • Vrfo + uo • V6) . 

ooi J 

( 12 ) 

Meanwhile, the magnetic moment fjL = and the 

energy E is defined as 

E = (13) 

The ion equilibrium distribution m can be rearranged 
as 


Fj = iV(V^)(^4||)-^e 


I 

—’7^- + 




(14) 


We focus on the ions perturbed distribution func¬ 
tion first. Using the properties of uq, one can show 
that Vuo = a;Ti?(ViiVC — VCVii) -b and 

uq ■ \/b = & • Vmo. Recalling the expression of $o in 
dH), we have 


VD-^i^ = ^{wl + n;i/2)V|| + 
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where a is short for T±/T\^, representing the anisotropy 
strength. It is also found 

(W||& -b Uo) ■ Vuo + Vuo • (W||&-b Uo) 

= (^-ba;Ti?^)va;T. (16) 


In view of the two equations above, the GK equation is 
reduced to 

. • , % E 'Z 

deShi — in'll sinOShi — i—rShi = ieJo — r-^:^S(j). (17) 

ujI ujI oE 

Here, loI = w\\/{qR) is the modified transit frequency 
and = kr5l, where 81 is the ions orbit width defined 

^«Zra[^f + M + ^^ + 2^lPT||M+^(l-a) . 
The disturbed distribution function can be easily solved 
as 


8E, = [I - Jl{krp^)] 
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In view of © and (fUll . one has 


OFn* 

OE 


J_ 9F^ _ _]_ pi 
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Generally for ions, the finite-Larmor-radius (FLR) ef¬ 
fect is taken into account by assuming ~ A <C 1 and 
the finite-orbit-width (FOW) effect is also considered by 
assuming ^ A. While for electrons, we have krPe — 0 
and n| ~ 0. The electron disturbed distribution function 
is simplified to 


= (19) 

® fe#0 

The dispersion relation of GAM is derived by us¬ 
ing quasi-neutrality condition, Sui = Sue, namely, 
/ d^wSEi = / d^wSEe- Inserting SEj into the quasi¬ 
neutrality condition gives birth to the dispersion relation. 
Before doing that, we need to simplify SFi by cutting 
off the coupling chains. Gonsidering S(j)n ^ {krP^Y8<f)Q, 
we take into account only the coupling between 8(j){j 
and by neglecting all high-order harmonics. Here, 

p® = Vt\\/^c ion Larmor radius. Consequently, we 

write 8Fi = 5F^ +8F^^e^''^. Eventually, using Eqs. (fT^ 
and m and the quasi-neutrality condition, S(j)k is ob¬ 
tained as [1 + t{1 + (^Z{C,))]5<j)k = —i^qkrp'^TkG{kM)S(j)o 
for k = ±1. For k = 0, we find 


S{M)S8^o - 

qkrP 


^5(-M)5</._i = 0. 
qkrP 

( 20 ) 


Finally, we obtain the following dispersion relation 
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m- 

(23) 


Here, C, is the normalized frequency defined as quiR/vT\\ 
and Z{C,) is the plasma dispersion function. Letting 
M = 0 and tr = 1 in the equation above reproduces 
the result in Ref. HI. For just cr = 1, the dispersion re¬ 
lation of GAM in a toroidally rotating isotropic plasma 
is recovered^. Assuming M = 0 and r = 0 in Eq. (EB) 
yields the result identical with the one in Ref. [13. 


IV. DISCUSSION AND CONCLUSION 

Explicit analytical solutions to the dispersion relation 
m for arbitrary C, are difficult to obtain. Here we are 
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restricted to the GAM with C ^ 1 in order to find an 
asymptotic solution, which requires large safety factor. 
Hence, analytical results below are expected to be accu¬ 
rate enough only when q is sufficiently high, i.e., q ^ 3 as 
shown in some previous literature^L^. Now we asymp¬ 
totically expand the plasma dispersion function = 

ia^/nexp (-C^)-C“Hl+C“^/2+3C“^/4-|-15C“®/8-|-) 

as usual and neglect all terms of order higher than 
and then we arrive at the following simplified 
dispersion relation: 


i ^ ^ = 0, (24) 


in which 

^ 2-|-4(7-|-(7^ (1 — (7)^ r 

Gl =- - \- T + ——— j —- -|- — 
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By assuming ( = q^tx + iq'yd, where fix and are 
both real and are the normalized frequency and damping 
rate of GAM, respectively, Eq. (IMl) yields the frequency 
of GAM as 


n 




(27) 


and the imaginary part of (l24l) gives the damping rate of 
GAM in a toroidally rotating anisotropic plasma 


Id = 




K 






(28) 


where y = 1 is adopted due to the factor that <C fix- 
The two equations above are the major results repre¬ 
senting the effects of toroidal rotation on the GAM in 
anisotropic plasmas. 


Eq. (l27l) illustrates the dependence of GAM frequency 
on the safety factor q, temperature ratio r, anisotropy 
O’, and the Mach number M. Eor general Mach number 
M ~ 0{1), by ignoring terms proportional to l/g^, GAM 
frequency is reduced to 


nl = 


+ T + 


4(1 -b r) 


-bM^ 




1 -b T 


(29) 


According to the frequency above, we can see that the 
coefficient of M'^ depends only on r and has no relation 
to anisotropy. The coefficient of relying on both cr 
and r is always positive and turns into 4 in the isotropic 
case. Thereby, the toroidal rotation always increases the 
GAM frequency^. 

Let us focus on the coupling effect of temperature ratio 
and the anisotropy by zeroing M. Eq. (1^ is reduced to 
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2 _ 
GO — 


1 

2 



(1^_ 

4(1 -b r) 


(30) 


In the isotropic limit with (7=1, the frequency above 
reproduces the classical one 7/4 -b r. It is found that 
when (7 < 3 -b 2 t, there is dGlQ^/dr > 0 and for a > 
3 -b 2r, dflQQ/dr is negative. That is to say, when a < 
3-b2r, the increased electron temperature (with fixed Ty) 
will increase the GAM frequency. When a is larger than 
3 -b 2r, which requires strong anisotropy since generally 
T > 1 in a realistic tokamak plasmas, the increasing of 
electron temperature decreases the GAM frequency. On 
the other hand, one can see (90c;o /da is always positive, 
meaning that the anisotropy always tends to enlarge the 
GAM frequency. More stronger anisotropy means more 
greater Sp± when py is fixed, leading to the increasing of 
the total perturbed pressure Sp\\ -b Sp±, which eventually 
determines the GAM frequency. Keeping terms to the 
order of 0(l/g^), Eq. (l27)l yields 


“G — “go 


1-b 


30 -b 16(7 -b (22 -b 8cr -b 2(7^)r -b 8 t^ -b 2(1 — a){T -b rtr — 2a) 
7^ (^2 -b 4(7 -b cr^ -b 4 t -b 


(31) 


This equation illustrates the approximate frequency of Eq. (1301) is not identical with the previous one reported 
GAM in anisotropic tokamak plasmas with arbitrary 
toroidal rotation in the large safety factor case. 
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in Ref. HO, which reads [see Eq. (23) in[ 
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Subtracted by the equation above, Eq. O leaves 

{a - 1)t(1 + O' + 2t) 


^Go “ ^Gi — 


4(1 + r) 


(32) 


(33) 


which indicates that our result becomes the same as the 
on in Ref. [13 only when r = 0 and/or cr = 1. Since 
Ref. HO used the standard bi-Maxwell distribution in 
which El and T|| were both assumed to be the func¬ 
tion of 'i/', dF^/d^ was neglected in the kinetic equation. 
However, a self-consistent distribution (THl) implies that 

= (1 — which can not be disregarded in the 

kinetic equation except for ct ~ 1. It is also noted that 
dF^/dfi only contributes to the first term on the right- 
hand side of Eq. (ITSl) . In the case of r = 0, current 
neutrality condition (J’') = 0 is used to obtain the dis¬ 
persion relation of GAM21 and only the second term on 
the right-hand side of Eq. (fT51) plays a role. Hence the 
dependence of Fq on the magnetic moment /r does not 
impact the final dispersion relation. 

The damping rate (l28ll is proportional to e~'^ and 
will be dramatically decreased by increasing Consid¬ 
ering that fix always increased by a and decreases with 
r when t < {a — 3)/2, we let cr = 3 and t = 0 to mini¬ 
mum which is about 6.75. As a result, the damping 
rate is nero zero and ignorable. For general anisotropy 
strength cr < 3, the temperature ratio r is usually greater 


than (cr — 3)/2. Then the GAM frequency increases as 
T grows up. It can be concluded that the collisionless 
damping rate is remarkably diminished by either r or cr. 

Geodesic acoustic mode (GAM) in a toroidally rotat¬ 
ing anisotropic tokamak is investigated by using gyro- 
kinetic (GK) equation. Self-consistent equilibrium dis¬ 
tribution function is derived according to the zeroth- 
order Fokker-Planck equation, leading to Tn = T||('0) 
and T± = — Bq{iIj)) for ions and Te = Te{ip) 

for electrons, where r ^ melmi is assumed. The per¬ 
turbed distribution function is analytically obtained from 
the GK equation with m = n = 0 for the perturbed 
electrostatic potential 5(f>. By focusing on the pass¬ 
ing ions and assuming small drift orbit radius, the gen¬ 
eral dispersion relation of GAM is derived. For large 
safety factor, the GAM frequency is given as Hqq = 
1/2 -I- cr-|- cr^/4-|- (1 — cr)^/[4(l -|- r)] to the leading order. 
The anisotropy strength a is shown to increase the GAM 
frequency and diminish the damping rate. The tempera¬ 
ture ratio increases the GAM frequency when it is larger 
than {a — 3)/2 and decreases the GAM frequency when 
T < (ct-3)/2. 
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